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Abstract 

In the framework of the left-right twin Higgs model, we study pair production of 
the neutral Higgs bosons at the International Linear Collider (ILC) and the CERN 

r- 1 \ 

LHC. We find that the production cross section of the process e + e~ — ► (f> h are 
at the level of several tens fb at the ILC, the production cross section of the (fPcfP 

> • 

pair and (fPh pair are at the level of several hundreds fb at the LHC . As long as 
' the neutral Higgs boson (fp is not too heavy, we conclude that its pair production 

. might be used to test for the left-right twin Higgs model at the LHC experiment 

O ■ 

' or in the future ILC experiment. 
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I. Introduction 



The Higgs mechanism is the heart of the standard model (SM) providing masses to 
gauge bosons via electroweak symmetry breaking (EWSB). However, the SM fails to 
explain the origin of the fermion mass and has naturalness problems. Many alternative 
new physics models with extended Higgs sectors are free from the above difficulties. The 
hunt for the Higgs bosons came to be one of the most important goals for present and 
future high energy collider experiments. Apart from the SM, neutral Higgs bosons appear 
in almost every scenario exploring new phenomena [1]. Pair production of neutral Higgs 
bosons at the CERN LHC, which will provide a way to test the Higgs boson self-coupling, 
may be sensitive to new physics [2, 3]. Many works have contributed to studies of the 
neutral Higgs pair production at the hadron collider in model independent [4], in SM 
[5-8] , and in new physics models beyond the SM, such as little Higgs models [9] , Randall- 
Sundrum-like models [10], top condensation models [11], supersymmetric models (SUSY) 
[12, 13] and models of universal extra dimensions (UED) [14]. 

The SM has been proved by all existing precise experimental data with its theoretical 
predictions beyond one-loop level being coincident with experimental observations. But 
in the SM the Higgs boson mass suffers from an instability under radiative corrections, 
which is called "hierarchy problem" [15]. Recently, the twin Higgs mechanism has been 
proposed as a solution to the little hierarchy problem. The Higgs bosons emerge as pseudo- 
Goldstone bosons once the global symmetry is spontaneously broken. Gauge and Yukawa 
interactions that break the global symmetry give masses to the Higgses. The twin Higgs 
mechanism can be implemented in left-right models with the additional discrete symmetry 
being identified with left-right symmetry [16, 17]. The left-right twin Higgs (LRTH) 
model contains the U(A)i x £7(4) 2 global symmetry as well as the gauged symmetry 
SU(2) L x SU(2) R x U(1)b-l- After Higgs obtained vacuum expectation values (/, /), 
the global symmetry £7(4)i x £7(4)2 breaks down to £7(3) i x £7(3)2, and the gauge group 
SU(2) R x U(1) B _ L breaks down to the SM £7(l) y . Thus, the LRTH model predicts 
the existence of the new particles, such as heavy gauge bosons, heavy scalars, and the 
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top partner T, which can generate rich phenomenology at present and in future collider 
experiments [17-21]. 

In the context of the LRTH model, pair production of the charged Higgs bosons (0 + , 
4>~) in the LRTH model at the ILC and LHC are studied in Ref.[21], but they did not 
consider production of the neutral Higgs bosons ((f) , h). As we know, so far production 
of the neutral Higgs pair at the LHC and the ILC in the LRTH model has not been 
considered, which is the main aim of this paper. 

Besides the SM-like Higgs boson h, there are two additional neutral Higgs bosons in 
the LRTH model, which are h% and 0°. The neutral Higgs boson h% is a possible dark 
matter candidate that only couples to the gauge bosons (including the SM gauge bosons 
7, Z, W , and the new gauge boson Z H ). The production cross section of h9, at the collider 
is very small and escapes the detector. Therefore, in this paper, we will not discuss the 
production of h9> at the ILC or LHC. The neutral Higgs boson 0° is a pseudoscalar that 
couples to both the SM fermions and gauge bosons. The neutral Higgs boson pair cf>°h can 
be produced via the processes e + e _ — > Z(Zh) — > <f>°h at the ILC, and via the partonic 
processes qq — > (f>°h(q = u,c,d, s,b), gg — > <p°h at the LHC, respectively. While the 
neutral Higgs pair 0°0° can only be produced via the partonic process gg — > and the 
t-channel partonic process 66 — > at the LHC. We calculate all above these processes. 
Our numerical results denote that, for m h = 120GeV, 120GeV < m^o < 180GeV and 
500GeV < f < 1500GeV: (i)the production cross section of <f)°h at the ILC with the 
center-of-mass (cm.) energy y/s = 500GeV is in the range of 0.92/6 ~ 20/6; (ii)the 
production cross section of (f)°h at the LHC with the cm. energy ^fs = lATeV is in 
the range of 34/6 — 306/6, and the main contribution comes from light quarks; (iii)the 
production cross section of at the LHC is in the range of 4/6— 122 fb, and the main 
contribution comes from the top quark loop. 

This paper is organized as follows. In Sec. II, we briefly review the essential features 
of the LRTH model. The relevant couplings of the neutral Higgs bosons to other particles 
and the feature of the decay for the neutral Higgs bosons 0° are also discussed in this 
section. In Sees. Ill and IV, we give our numerical results for pair production of neutral 
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Higgs bosons predicted by the LRTH model at the ILC and LHC, respectively. Our 
conclusions are given in Sec. V. 

II. The LRTH Model 

The LRTH model was first proposed in Ref.[16] and the details of the model as well 
as the particle spectrum, Feynman rules, and some phenomenology analysis have been 
studied in Ref.[17]. Here we will briefly review the essential features of the model and 
focus our attention on the neutral Higgs bosons. 

The LRTH model is based on the global U(A)i x U(A) 2 symmetry with a locally 
gauged subgroup SU(2) L x SU{2) R x U(1) B -l- Two Higgs fields, H = {H L ,H R ) and 
H = (H L ,H R ), are introduced and each transforms as (4,1) and (1,4), respectively, 
under the global symmetry. H^ R (Hl )R ) are two component objects which are charged 
under SU(2) L and SU(2) R , respectively. For the gauge couplings g 2 L and g 2R of SU(2) L 
and SU(2) R , the left-right symmetry implies that g 2 L = 92R = 92- 

The U(A)i [C/(4) 2 ] group is spontaneously broken down to its subgroup U(3)i [f7(3) 2 ] 
with nonzero vacuum expectation value (VEV) < H > = (0, 0, 0, /) [< H > = (0, 0, 
0, /)]. The Higgs VEVs also break SU(2) R x U(1) B -l down to the SM U(1) Y - After 
spontaneous global symmetry breaking by / and /, three Goldstone bosons are eaten by 
the new gauge bosons and Z H . After the SM electroweak symmetry breaking, the 
three additional Goldstone bosons are eaten by the SM gauge bosons W ± and Z. 

The fermion sector of the LRTH model is similar to that of the SM, with the right- 
handed quarks (u R , d R ) and leptons (l R , v R ) form fundamental representations of SU {2) R . 
In order to give the top-quark mass of the order of the electroweak scale, a pair of vectorlike 
quarks Ql and Q R are introduced. The mass eigenstates, which contain one the SM top 
quark t and a heavy top partner T, are mixtures of the gauge eigenstates. Their masses 
are given by 

m 2 t = 1 -{M 2 + y 2 f 2 - N t ), M 2 = ^(M 2 + y 2 f 2 + N t ). (1) 

where N t = ^J(y 2 f 2 + M 2 ) 2 - y 4 f 4 sin 2 2x with x = v/V2f, in which v = 2AQGeV is the 
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scale of the EWSB. Provided M T < f and that the parameter y is of order one, the top 
Yukawa coupling will also be of order one. The parameter M is essential to the mixing 
between the SM top quark and its partner T. 

According the symmetry-breaking pattern discussed above, with certain reparametriza- 
tions of the fields, there are left with four Higgs bosons in the LRTH spectrum that couple 
to both the fermion sector and the gauge boson sector. They are one neutral Higgs bosons 
0°, a pair of charged Higgs bosons ± , and the SM-like physical Higgs h. In addition, 
there is an SU(2) L doublet h = (hi^h^) that couples to the gauge boson sector only 
(including the SM gauge bosons 7, Z, W , and the new gauge boson Z H ). The lightest 
particle in h, typically one of the neutral components, is stable, and therefore constitutes 
a good dark matter candidate. 

These neutral Higgs bosons can couple to each others, and also can couple to the 
ordinary fermions, ordinary gauge bosons, new top quark T, and new gauge boson Zh- 
The couplings expression forms which are related our calculation, are shown as [17] 



Where pi, p 2 , and p 3 refer to the incoming momentum of the first, second and third par- 
ticles, respectively. Uj and di represent the upper- and down- type fermions, respectively. 
Sw = sinOw, Cw = cos9w, and 9w is the Weinberg angle. At the leading order of 1/ /, 
the sine values of the mixing angles and can be written as 
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Cl and Cr are the cosine values of the mixing angles «i and or, respectively. Pl(r) = 
(1 =f 75)/2 is the left (right)-handed projection operator. 

In the framework of the LRTH model, the mass of the neutral Higgs boson 0° can 
be anything below / here we consider another possibility, in which the mass is around 
150GeV [17]. Similar to the SM Higgs boson, 0° can decay to 77 through the top-quark 
loop and heavy top-quark loop. But unlike the SM Higgs boson, in the LRTH model, 
the light neutral Higgs boson 0° is a pseudoscalar boson, due to its pseudoscalar nature, 
there is no 4>°WW and 4>°ZZ couplings at tree level. So, the one-loop SM gauge boson 
contribution to 0°77 is zero. In general, the light neutral Higgs boson 0° decays into 66, 
cc, t + t~, gg and 77. Now we discuss the branching ratios for the possible decay modes 
of 0°. The decay width of 0° — > // is proportional to the square of the corresponding 
Yukawa coupling, with an additional suppression factor of z/ 2 /(2/ 2 ) comparing to that 
of the SM Higgs boson. The concrete expressions of the decay widths for the different 
decay channels are given as follows: 

r(0 u ^66) = g ^ (i-HK°) 2 . 

3G F m 0O z/ 2 m 2 
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r(0° - r+o = 



Where m b , m c , and m T are the masses of the SM fermions 6, c and r, respectively. The 
index / corresponds to q and I (q = quark, I = lepton). Where N[ = 1,3 for / = l,q, 
respectively. Qf is the charge of the fermion /. Similar with Ref.[22], the function 
can be written as: 

Af = 2r f [l + (1-T f )f(r f )}. (5) 

where 77 = 4m 2 /m|o. In Ref. [22], the function /(t/) has two parts corresponding to 
the Tf > 1 and 7/ < 1 two conditions. In our numerical estimation, we have neglected 



the contributions of the light fermions. Therefore, in the LRTH model, there is r^m 

t(T) / m 0° 



4mLi/m^ > 1, and the function /(r/) is given by 



fir) = arcsin 2 -j=. (6) 
where , r g , and /(r g ) in Eq.(4) are defined the same as , 77 and /(r/), but only for 
quarks. 

Using above partial widths of the neutral Higgs boson 0°, its total width Y can be 
approximately written as 

r = ^bi + r c g + r r + T - + r 77 + r gg . (7) 
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Figure 1: The branching ratios of the neutral Higgs boson </>° for different decay modes 
as functions of the free parameter / for M = 150GeV, m^o = 120GeV. In 



order to see the trend clearly, we have multiplied Br ((f) — > cc), Br( 



r 1 r 



and Br I 



77) by the factors 10, 20, and 300, respectively. 



We summed up our numerical results of the branching ratios of the neutral Higgs 
boson 0° for different decay modes Br ((ft ) in Fig.l. To get the numerical results, the 
SM parameters involved are taken as mj, = 4.8GeV, m c = 1.25GeV and m T = 1.78GeV 
[23]. In Fig.l, we plot Br((f>°) as a function of free parameter / for M = 150GeV 



and m^o = 120GeV. One can see from Fig.l that the decay branching ratios of (f>° 
are sensitive to the parameter /. If we assume that the parameter / is in the range of 
500GeV ~ lSOOGeV, the value of the branching ratio Br((f>° — > bb) is in the range of 
14% — 55%, and the branching ratio Br((f>° — > gg) is in the range of 38% — 85%. The 
values of Br((f>° — * cc), Br((j) — > t + t~) and Br(<f)° — > 77) are much smaller than those 
of Br (<f) — > 66) and Br ((f) — > gg). Therefore, in order to see the trend clearly, in Fig.l 
we have multiplied them by 10, 20 and 300, respectively The real numerical results 
are 5r(0° -> cc)= 0.9% - 3.8%, Br((f)° -> r+T")= 0.6% - 2.6%, and 5r(0° -> 77) = 
0.09% — 0.2%. Our numerical results agree quite well with Ref.[17], in that the branching 
ratio Br ((f) — > 77) is roughly same as _Br(/i — > 77) for rrih = m^o. 

III. Pair production of neutral Higgs bosons at the ILC 

In many cases, the ILC can significantly improve the LHC measurements. If a Higgs 
boson is discovered, it will be crucial to determine its couplings with high accuracy, to 
understand the so-called mechanism of EWSB [24]. The high resolution profile determi- 
nation of a light Higgs boson (mass, couplings, self-couplings, etc.) can be carried out 
at the ILC, where clear signals of Higgs events are expected with backgrounds that can 
be reduced to a magnitude level. With the LHC guidance, the ILC, which is currently 
being designed, will further improve our knowledge of the Higgs sector if that is how 
nature decided to create mass [24]. It was demonstrated in Ref.[25] that physics at the 
LHC and at the ILC will be complementary to each other in many respects. So far, 
many works have been contributed to studies of the neutral Higgs boson pair production 
at the ILC, in the SM [26-28] and in new physics beyond the SM [29-32]. 

From the discussions given in Sec. II, we can see that the neutral Higgs boson pair 
cannot be produced exclusively at the ILC because 0°0°can not couple with gauge 
boson Z or Zh- However, the neutral Higgs boson pair (f)°h can be produced via the 
processes e + e~ — > Z(Z H ) — > <f)°h at the ILC. The Feynman diagrams of the process 
e + (pi)e~ (p 2 ) — > <t>°{pz)h(pi) are shown in Fig.l. 
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Figure 2: Feynman diagrams for the process e + e~ — ► <fr°h. 
At the leading order, the production amplitude of the process can be written as 

M- M z I M Zh 
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where P12 is the momentum of the propagator, which is the sum of the incoming momen- 
tums pi and p 2 . With the above production amplitudes, we can obtain the production 
cross section directly. 

From the above discussions, we can see that, except for the SM input parameters 
a = 1/128.8, S w = V0.2315, m z = 91.1876GeV, m h = 120GeV [23], the cross section a 
of pair production for the neutral Higgs boson <ft°h at the ILC is dependent on the model 
dependent parameters / and m^o. In our numerical estimation, we will assume that the 
values of the free parameters / and m^o are in the ranges of 500GeV — lSOOGeV^ and 
lQQGeV - 180GeV, respectively. 

In Fig. 3, we plot the production cross section a of the process e + e~ — *> <f>°h as a 
function of the scale parameter / for the cm. energy y/s = 500Ge^, rrih = 120GeV and 
three values of m^o . We can see that a is sensitive to the scale parameter / and the mass 
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parameter m^o. For 500GeV < f < 1500GeV and UOGeV < m o < 180GeV, its value is 
in the range of 0.92/6 — 20/6. According to an update of pa—fo^^at 2006 [Sag 
one can see that, an integrated luminosity of 500/6 -1 should be achieved in the first four 
years of running after one year of commissioning. Therefore, if we assume the integrated 
luminosity for the ILC is 500/6 -1 , there will be 10 2 — 10 4 4>°h events to be generated at 
the ILC. 




—\ 1 1 > 1 > 1 > 1 — 
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Figure 3: The production cross section o of e + e~ — > <f)°h as a function of the parameter 
/ for three values of m^o, nih = 120Ge^, and the cm. energy y/s = 500GeV^. 

From the discussions given in Sec. II, we can see that the possible decay modes of the 
neutral Higgs boson 0° are bb, cc, t + t~, gg and 77. The SM-like neutral Higgs boson 
h has similar decay features with those of 0°. Therefore, the signatures of neutral Higgs 
boson pair <p°h is similar to those of the neutral Higgs boson pair at the high energy 
colliders. From the numerical results given in Sec. II, one can see that, for the masses 
m^o < 180GeV, the possible signals of (j)°h can be seen as four b quarks, 

e + e- -> (jPh -> bibb. (9) 

The production rate of the bbbb final state in the LRTH model can be easily estimated 
1 Thanks to the referees for offering this reference to us. 
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using the formula a s = a x Br((p° — > bb) x Br{h — > 66). If we assume the integrated 
luminosity £ int = 500/6 -1 for the ILC with the cm. energy y/s = bWGeV , then there 
will be 9 — 3.0 x 10 3 6666 events to be generated at the ILC, which is significantly larger 
than that for the SM Higgs boson pair production process e + e~ — > hh — > 6666 [26-28]. 
Therefore, we hope that by using very efficient yU-vertex detectors to tag the 6 quark jets, 
we might detect the possible signatures of the neutral Higgs boson 0° via the process 
e + e~ — > <p°h in the future ILC experiments. Certainly, detailed confirmation of the 
observability of the signals generated by the process e + e~ — > Z(Z H ) — > would require 
Monte-Carlo simulations of the signals and backgrounds, which is beyond the scope of 
this paper. 

IV. Pair production of neutral Higgs bosons at the LHC 

The LHC has a good potential for discovery of a neutral Higgs boson. Now we look 
at pair production of the neutral Higgs bosons predicted by the LRTH model at the 
LHC. From the above discussions, we can see that both the (/>°0° pair and <p°h pair can 
be produced at the LHC . In this section, we will consider both of these cases. 

A. pair production 

First, we study production of the neutral Higgs boson pair at the LHC . At 
the LHC , the neutral Higgs boson pair can be produced through two mechanisms. 
One is loop-induced production via gluon fusion (gg — > and the other is from the 

t-channel quark- ant iquark annihilation (qq — > The relevant Feynman diagrams 

are shown in Fig A. Considering the couplings of the neutral Higgs boson 0° to the SM 
fermions are proportional to the factor of m q j f and the smallness masses of the quarks 
q — u,c, d, and s, we have neglected their contributions to production of the neutral Higgs 
boson pair 0°0°. 

In this paper, we calculate all production channels for the neutral Higgs boson pair 
<f)°(f) at the LHC, as shown in Fig A, including triangle diagrams, box diagrams and tree- 
level diagram. Each loop diagram is composed of some scalar loop functions, which are 
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calculated by using LoopTools [34]. The hadronic cross section at the LHC is obtained by 
convoluting the partonic cross sections with the parton distribution functions (PDFs). In 
our numerical calculation, we will use CTEQ6L PDFs for the gluon and quark PDFs [35] . 
The renormalization scale hr and the factorization scale \xp are chosen to be \xr — fip — 
2m0o. Because the calculation of the loop diagrams are too tedious and the analytical 
expression are lengthy, we will not present those here. 




b 1 

(c) 

Figure 4: One-loop Feynman diagrams for the subprocess gg — > 0°0°(a,b) and tree- 
level Feynman diagram for the subprocess bb — ► 0°0°(c) in the LRTH 
model. The diagrams obtained by exchanging the two gluons or exchanging 
the two Higgs bosons are not shown here. 

It is obvious that the production cross section a of the neutral Higgs boson pair (f)°(f) 
at the LHC are dependent on the model dependent parameters /, m^o, and M. Similar 
to the calculation at the ILC, we assume that the values of the free parameters / and m^o 
are in the ranges of SOOGeV^ — lSOOGeV^ and lOOGeV — 180GeV, respectively. Besides, 
we assume the mixing parameter M is in the range of lOOGeV — 200GeV. Our numerical 
results are summarized in Figs.5 and 6. 

To see contributions of the different partonic processes to the total hadronic cross 
section, we plot the total and partial hadronic cross sections for different partonic processes 
as functions of the scale parameter / for the parameters M = lOOGeV and m^o = 120GeV 
in Fig.5. We see from Fig.5 that production of the neutral Higgs boson pair is 
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dominated by the partonic process gg — > induced by the top-quark loop diagrams. 
For M = lOOGe^, m o = UOGeV, and 500GeV < f < 150GeV, the value of the total 
production cross section is in the range of 4/6 ~ 122/6, and the value of the production 
cross section coming from the top-quark loop diagrams is in the range of 1.5/6 — 105/6. 
This is because the contributions of the box diagrams are generally much smaller than 
those of the triangle diagrams, and furthermore the coupling hti is much larger than the 
coupling hTT or the coupling 0°66. If we assume the integrated luminosity £i nt = 100/6" 1 
for the LHC with the cm. energy yfs = 14TeV, then there will be 4 x 10 2 — 1.22 x 10 4 
events to be generated at the LHC. 



140 




600 800 1000 1200 

f(GeV) 



Figure 5: The total and partial hadronic cross sections for different partonic processes 
as functions of the free parameter / for the parameters M = lOOGeV and 
m^o = UOGeV. 

In order to see the effects of the mass parameter m^o on the total cross section a, we 
plot a as a function of m^o for / = 500GeV and three values of the mixing parameter M 
in Fig.6. One can see from Fig.6 that the total cross section a is sensitive to the mass 
parameter m^o, while is not sensitive to the mixing parameter M. This is because M 
is introduced to generate the mass mixing term Mq^qn, which is included in the gauge 
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invariant top Yukawa terms allowed by gauge invariance. From the relevant Feynman 
rules we can see that, the mixing parameter M does not influence the production cross 
section a of the neutral Higgs boson 0° too much. For / = 500Ge^, M = 200GeV, and 
m^o = lOOGeV — 180Ge^, the total cross section a is in the range of 16/6 — 253/6. 
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Figure 6: The total production cross section a as a function of free parameter m^o for 
three values of mixing parameter M. 

B. (fPh pair production 

Now we consider production of the neutral Higgs boson pair <j)°h at the LHC. At the 
LHC, the neutral Higgs boson pair <p°h can be mainly produced through two mechanisms: 
(i) qq — ► (j)°h, where q — u,d, c, s, b; (ii) the loop-induced gluon fusion process gg — > <j>°h. 
The relevant Feynman diagrams are shown in Fig. 7. 

Using the relevant Feynman rules, we can write the invariant amplitude for the par- 
tonic process q{j>i)q{j>2) ( P (P3)h( y p4) as 

M 2 (q) = M 21 (q), for q = u, c 

M 2 (q) = M 22 (q), forq = d,s 

M 2 {q) = M 22 {q) + M 23 (q), forq = b (10) 
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Figure 7: Tree-level Feynman diagrams for the process qq — > <f)°h(q = u, d, c, s, b) (a,b) 
and one-loop Feynman diagrams for gg — > (f)°h (c,d,e) in the LRTH model. 

For the s-channel partonic processes qq — > Z(Zu) — > ((Ph (q = u and c), the invariant 
amplitude can be written 
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For the s-channel partonic processes qq — > Z(Zh) —* <jPh (q = d, s and 6), the invariant 
amplitude can be written 
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For the t-channel partonic process bb — > <p°h as shown in Fig. 7b, the invariant ampli- 
tude can be written 



M 23 (g) = —nr-v(p 2 )— 2 2 l^u{pi). 



V2vf p\z - ml 

Where p 13 = p\ — p 3 . Considering the couplings of the neutral Higgs boson 0° to the SM 
fermions are proportional to the factor of m q /f and the smallness masses of the quark 
q — u,c, d, and s, we have neglected their contributions to production cross section of 
the neutral Higgs boson pair <p°h via the t-channel process in our calculations. When 
we calculate the loop diagrams Fi gs. 7 (c)-7(d), and Fig.7e, we will use the same method 
with FigsA(a) and 4(b). 

To see contributions of the different partonic processes to the total hadronic cross 
section, we plot the total and partial hadronic cross sections for different partonic processes 
as functions of the parameter / for m^o = m h = l20GeV and M = 150GeV in Fig.8. We 
see that the production cross sections of the neutral Higgs bosons (f)°h mainly come from 
the contributions of the light quarks (u, d, c, s) through the s-channel Z exchange and Zh 
exchange. Our numerical results show that, the contributions coming from the partonic 
processes gg — > <jPh [including Figs.7(c)-7(e)} to total production cross section are at the 
orders of 10~ 5 /6 — 10~ 1 fb, which are much smaller than those of the tree-level processes. 
This is because the Yukawa couplings depend sensitively on the free parameters M and 
/. The parameter M is very smaller than the scale parameter /. So, although the gluon 
fusion get an enhancement due to large parton distribution functions, the contribution 
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of the gluon fusion process is suppressed by the order of [Mj f) 4 [21]. Thus, in Fig. 8, 
we did not show the line corresponding to the value of the production cross section 
contributed by the gg fusion. The value of the production cross section of the neutral 
Higgs bosons (j)°h is insensitive to the mixing parameter M. For m^o = = 120GeV 
and 500GeV < / < 150GeV, its value is in the range of 34/6 — 306/6, the partial value of 
the total production cross section coming from light quarks contributions is in the range 
of 31/6 — 281/6. If we assume the integrated luminosity £i n t = 100/6 -1 for the LHC 
with the cm. energy yfs = 14TeV, then there will be 3.4 x 10 3 — 3.1 x 10 4 <p°h events 
generated at the LHC. 




600 800 1000 1200 1400 



f(GeV) 

Figure 8: The total and partial hadronic cross sections for different partonic processes 
as function of the parameter / for m^o = m/ t = 120GeV and M = 150GeV. 

Similar to those of the discussions for neutral Higgs boson pair (/>°</>° production, we 
plot a as a function of free parameter / for = 120GeV, M = 150GeV and three values 
of m^o in Fig.9. One can see from Fig.9 that the total cross section a is sensitive to mass 
parameter m^o. For / = 500GeV and 120GeV < m^o < 180GeV, its value is in the range 
of 101/6- 306/6. 

From the above discussions, we can see that the decay features of (ft are similar to 
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those of the SM-like neutral Higgs boson h, as far as decays into bb and 77 are concerned. 
Therefore, when we analyze the signatures of the neutral Higgs boson pairs from the 
LRTH model at the colliders, we will take the 0°0° pair, for example. 

300- 
250- 
200- 

n 

^0 150- 
100 - 
50- 

600 800 1000 1200 1400 

f(GeV) 

Figure 9: The total production cross section as a function of free parameter / for = 
UOGeV, M = 150GeV and three values of m o. 

In most of the parameter space of the LRTH model, the main decay modes of 0° 
are gg and bb. However, the final states gggg and bbbb induced by pair production of 
the neutral Higgs boson 0° at the LHC have large QCD backgrounds and thus are 
insignificant for 0° discovery. If we assume that one of the neutral Higgs boson 0° decays 
to bb and the other decays to 77, then pair production of the neutral Higgs boson 0° at 
the LHC can give rise to the 6677 final state, and the production rate of the 6677 final 
state can be easily estimated using the formula a s = o x Br(<f)° — ► bb) x Br '(0° — > 77). If 
we assume the integrated luminosity £i nt = 100/6 -1 for the LHC with the cm. energy 
^/s = 14TeV, then there will be several hundreds of 6677 events to be generated at the 
LHC. Furthermore, the narrow 77 peak can be reconstructed to distinguish the signal 
from the backgrounds. Detailed analysis of the signals and the relevant backgrounds 
about this kind of the final state has been given in Ref . [36] . 
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V. Conclusions 



The twin Higgs mechanism provides an alternative method to solve the little hierarchy 
problem. The LRTH model is a concrete realization of the twin Higgs mechanism. In 
this paper, we discuss the possible decay modes of the neutral Higgs boson 0° predicted 
by the LRTH model and consider its pair production at the ILC and LHC via suitable 
mechanisms. 

At the ILC, we study production of the neutral Higgs boson pair <p°h via the processes 
e + e~ — > Z(Zh) — > <f>°h. Our numerical results show that, for m^o = m^ = 120GeV and 
500GeV < f < 1500GeV, the total production cross section of neutral Higgs boson pair 
(f)°h at ILC is in the range of 0.92/6 — 20/6. If we assume the integrated luminosity 
£int = 500/6~ 4 for the ILC with the cm. energy yfs = 500GeV, there will be 10 2 — 10 4 
<f)°h events to be generated at the ILC. If we assume that the neutral Higgs bosons 0° 
and h both decay to bb, then the process e + e~ — > (j)°h can give rise to the 6666 final state. 
There will be 9 — 3.0 x 10 3 bbbb events to be generated at the ILC. Owing to the bbbb 
events, we might detect the possible signatures of the neutral Higgs boson 0° via the 
processes e + e~ — > Z(Zh) — > (f>°h in the future ILC experiments. 

At the LHC, we study production of the neutral Higgs boson pairs 0°0° and <jPh. First, 
we study production of the neutral Higgs boson pair 0°0° via the processes gg — > 0°0° 
and qq — > 0°0°. Our numerical results show that, for M = lWGeV, m^o = 120GeV 
and 500GeV < f < 1500GeV, the value of the hadronic cross section a^o^o is in the 
range of 4/6 — 122/6, which mainly comes from the contributions of the top-quark loop. 
Then we study production of the neutral Higgs boson pair (f>°h via the processes qq — > 
<p°h(q = u, c, d, s, 6) and gg — > <f>°h. Our numerical results show that, for M = 150GeV, 
m^o = m h = 120GeV and 500GeV < f < 1500GeV, the value of is in the range of 
34/6 — 306/6, of which about 91% of the contributions comes from light quarks u, d, c, s. 
If we assume the integrated luminosity £i n t = 100/6 -1 for the LHC with the cm. energy 
y/s = lATeV, then there will be 3.4 x 10 3 — 3.1 x 10 4 (jPh events to be generated at the 
LHC. If we assume that one of the neutral Higgs bosons 0° and h decays to 66 and the 
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other decays to 77, then the processes pp — > + X and pp — > <j)°h + X all can give rise 
to the 6677 final state. There will be several hundreds and up to thousands of 6677 events 
to be generated at the LHC with the cm. energy y/s = lATeV and £ int = lOO/fr -1 . 
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